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Abstract. We discuss the band-gap structure and the integrated density of 
states for periodic elliptic operators in the Hilbert space L 2 (S. rn ), for m > 2. 
We specifically consider situations where high contrast in the coefficients leads 
to weak coupling between the period cells. Weak coupling of periodic systems 
frequently produces spectral gaps or spectral concentration. 

Our examples include Schrodinger operators, elliptic operators in divergence 
form, Laplace-Beltrami-operators, Schrodinger and Pauli operators with peri- 
odic magnetic fields. There are corresponding applications in heat and wave 
propagation, quantum mechanics, and photonic crystals. 



Introduction 

We consider periodic elliptic partial differential operators of the second order 
with high contrast in (some of) the coefficients. The classes of operators discussed 
here include Schrodinger operators, magnetic Schrodinger and Pauli operators, 
the acoustic operator, Laplace-Beltrami operators, and Maxwell operators. We 
are mainly interested in the question whether high contrast in the coefficients 
leads to weak coupling between the period cells and to spectral concentration 
and/or spectral gaps. For most of the paper, we assume a simple geometry like 
the one shown in Figure 1, for which the strongest results are obtained. We 
include some remarks on more general situations, like sponges in dimensions 3 
and higher. Our main theme is to show that simple topological assumptions lead 
to spectral consequences for a wide class of — rather different — elliptic operators. 
We do not discuss the question of absolute continuity here. 

The paper is organized as follows. Section 1 contains some basic definitions, 
assumptions and background results concerning monotone convergence of qua- 
dratic forms. In Section 2, we briefly report on some "classical" results on the 
convergence of Schrodinger operators with high barriers. Similar results for strong 
magnetic barriers in the Schrodinger and Pauli case are discussed in Section 3. 
Section 4 is devoted to divergence type operators —V ■ (1 + Axrz)V, where Q is 
a periodic, open subset of R m , with m > 2. There is some superficial similarity 
between divergence type operators and Laplace-Beltrami operators. The exis- 
tence of spectral gaps for periodic Laplace-Beltrami operators with a conformal 
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Figure 1 . The fundamental period cell Q and the set M = Q n M. 

factor in the metric tensor is discussed in Section 5. We also refer to construc- 
tions of a purely "geometric" nature that yield Laplace-Beltrami-operators with 
gaps. We conclude in Section 6 with some references to the work of A. Figotin, 
P. Kuchment and others on the Maxwell operator and applications to photonic 
crystals. 

1. Preliminaries and basic assumptions 

Throughout the present paper, we will always assume that f2 C W n is open and 
periodic w.r.t. some discrete lattice T; for most purposes it is enough to consider 
P = Z m . Q denotes a fundamental period cell of T; we may take Q = (0, l] m 
if T = Z m . The complement of Q is denoted as M, a closed and periodic set. 
We will study two-component media with material parameter A = 1 on M, and 
A > 1 on fi. The strongest results are obtained in the "model situation" of 
Figure 1 where Q contains the boundary of the period cell Q. 

Our basic Hilbert spaces are TC = L 2 (R m ) and 7i 1 (M m ), the Sobolev space of 
L 2 -functions that have weak derivatives in L 2 . 

All our operators are semi-bounded, self-adjoint operators T acting in some 
L 2 -space; they are obtained from some closed, semi-bounded quadratic form t 
via the usual representation theorem ([K, RS-I]) and satisfy D(T) C D(t) and 
the relation 

t[ u , v] = {Tu, v) , u e D(T), veD(t). (1.1) 

As usual, we say T 2 > 7\ iff the associated quadratic forms satisfy -D(t 2 ) C D(ti) 
and t 2 [u, u] > t±[u, u], for all u e D(t 2 ). We also write t[u] := t[u, u] in the sequel. 
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For a sequence of non-negative, self-adjoint operators T n in Ti we say that T n 
converges to a limit Too in the strong resolvent sense iff 

(T n + 7)- 1 /^(T 0O + 7)- 1 /, feH; (1.2) 

for convergence in the norm resolvent sense we require uniform convergence for 
all / e W s.th. \\f\\ < 1. 

Most of our results hinge upon the following convergence theorem for quadratic 
forms ([RS-I], [S]): 

1.1. Theorem. Let (T n ) be a monotonic sequence of non-negative, self-adjoint 
operators. Then there exists a self-adjoint operator > s.th. T n — > Too in 
strong resolvent sense. 

For an increasing sequence of quadratic forms < U < t„ + i with domains Q n 
the limiting form domain is given by Qoo = {uG H^^Qn; sup t n [u] < oo}. Here 
it may happen that Too acts in a strictly smaller Hilbert space than the T n . In 
this case, (T^ + J) -1 has to be complemented by the zero operator on a suitable 
subspace of n. For a decreasing sequence it is shown in [RS-I, S] that there exists 
a largest closable form that is smaller than t„, for all n. 

In a periodic situation, strong resolvent convergence often implies convergence 
of spectral densities. For the definition of the integrated density of states (i.d.s.) 
we refer to [RS-IV], [PF] or [DIT]. In our context, the i.d.s. is a monotonically 
increasing function F : R — > R that tends to zero at — oo. Such functions have 
at most a countable number of discontinuities. The associated Borel measure \x 
is then called the density of states measure. We will say that a sequence (//„) of 
density of states measures converges weakly to fioo if F n (t) — > F^t) at all points 
of continuity t of F^. The proof of the subsequent proposition is elementary if 
one uses the definition of [RS-IV] for the i.d.s. 

1.2. Proposition. In addition to the assumptions of Theorem 1.1, suppose that 
the operators T n are periodic ( with respect to the the same lattice V ), and that a 
(finitely-valued) i.d.s. F n exist for T n as well as for Too. Then the corresponding 
measures fx n converge weakly to fioo- 

It is well-known that norm resolvent convergence implies spectral convergence 
(cf. [RS-I; Ch. VIII]) on any compact interval of the real line. In particular, we 
have the following simple consequence for the existence of spectral gaps: suppose 
T n — > Too in norm-resolvent sense. If Too has a gap around E e R, then the 
T n will also have a gap around E, for large n. More precisely, suppose that 
(a, b) H a (Too) = and let e > 0. Then (a + e, b - e) n o~(T n ) = 0, for all n > n e . 

2. SCHRODINGER OPERATORS 

As a warm-up, we illustrate the above convergence theorems by some simple re- 
sults on Schrodinger operators of the form — A + XV (x), where V is non-negative. 
It has been a "folk theorem" for quite some time that —A + Xxb converges in 
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strong resolvent sense to the Dirichlet Laplacian on the complement M m \ B, 
under suitable assumptions on B. 

Suppose V(x) > is a periodic, real-valued (continuous or measurable and 
bounded) potential s.th. {i G M m V(x) ^ 0} = Q. We let 

H x = -A + \V(x), A>0, (2.1) 

with associated quadratic forms f)\[u] = ||V-u|| 2 + A / V^w) 2 , for u G H}. 

To describe the limit, we will need here (and in the following sections) the 
Dirichlet Laplacian on the closed set M = W n \Vl ) denoted as — Am ([Hed, HZ]). 
The operator — A M is constructed from the quadratic form ||Vm|| 2 with domain 

Hl(M) := {u G H^R" 1 ) | u(x) = a.e. in M c }. (2.3) 

For sufficiently regular boundary of M, —Am agrees with the usual Dirichlet 
Laplacian. In the general case, it is shown in [HZ] that there exists a Borel 
set M*, M° C M* C M, such that —Am acts as a self- adjoint operator in the 
Hilbert-space L 2 (M*). It is now easy to see that 

{u G H 1 ^™) | sup f)>] <oo} = nl(M), (2.4) 

A>0 

and so Theorem 1.1 implies that H\ converges in strong resolvent sense to —Am- 
(To be more precise, the resolvent of — Am acts in L 2 (M*) and has to be com- 
plemented by the zero operator on L2(M m \ M*) in Eqn. (1.1).) Up to this point, 
periodicity was not required. If we now assume that Q is periodic, then for all 
operators involved the i.d.s. is well-defined, and we may also conclude that the 
i.d.s. converges. 

As was observed in [HH-1], one can actually do better and upgrade to norm 
resolvent convergence. The proof combines uniform exponential decay estimates 
for the resolvent kernels with Schur's Lemma and local compactness of —A. 

The above convergence results do not require Q to be connected. In our "model 
case" of Figure 1 where Q contains a neighborhood of dQ, —Am is a countable 
direct sum of copies of — Amhq- Therefore, the spectrum of —Am is a discrete 
set with each point in the spectrum an eigenvalue of infinite multiplicity. We 
therefore see that the spectrum of H\ concentrates at a discrete set of points, 
and an arbitrarily large number of spectral gaps opens up, as A — > oo. 

3. Magnetic Schrodinger and Pauli operators 

Here we briefly mention some results of [HH-1, 2], [HN], [B-1,2] on strongly 
coupled periodic magnetic fields. 

(1) Let a G C 1 (M m ,lR m ) be a magnetic vector potential and B = da the 
associated magnetic field. We assume that B is periodic w.r.t. the lattice T, and 
that B(x) = for all x G M, while B(x) ^ for all x G Vt. Note that periodicity 
of a is not required. The magnetic Schrodinger operator is then given as [CFKS] 

H(\a) = (-iV - Aa) 2 , (3.1) 
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and obtained via the closure of the forms 

(H(Xa)u,u) = ||(-iV - \a{x))uf, u G CfiW 1 ). 

Evidently, there is no monotonicity w.r.t. A. However, in any open subset U 
of Q where a fixed entry Bij(x) of the magnetic tensor is bounded away from 
zero, the Avron-Herbst-Simon bound [AHS] gives a local lower bound for the 
quadratic form, expressing the barrier-like effect of a strong magnetic field. If 
the vector potential a vanishes on M, the above results for —A + nxui I 1 ~^ °°, 
can be combined with the Feynman-Kac-Ito-formula to obtain strong resolvent 
convergence of H(Xa) to —Am, as A — > oo ([HH-1]). For given magnetic field B, 
it depends on the topology of the set Q and on flux conditions whether a gauge 
exists s.th. da = B and a(x) = a.e. on M. 

Under the above conditions, we also obtain an upgrade to norm-resolvent con- 
vergence, as before, and get the same spectral consequences as for —A + Xxn- 

(2) The situation where a(x) is allowed to have non-zero values on M has been 
analyzed in detail in [HN], using some basic homology theory and de Rham's 
Theorem. It can be shown that, under suitable assumptions, the norm differ- 
ence of the resolvents of H(Xa) and Hm{Xo) tends to zero, as A — > oo, where 
Hm(Xo) = (— iV — Aa) 2 , acting in /^(M*) with Dirichlet boundary conditions. 
As a consequence, one can conclude that the spectrum of H(Xa) approaches a 
periodic or a quasi-periodic function of A, as A — > oo. The precise result involves 
the flux of B through the connected components of (M°) c , the complement of 
the closure of the interior of M. 

(3) For the Pauli operator in M 2 , we have the pair of operators 

H±(Xa) = H(Xa) t B, (3.2) 

both acting in L 2 (R 2 ). Here supersymmetry (cf. [CFKS]) implies 

a(H + (Xa))\{0} = o-(H_(Xa))\{0}. (3.3) 

Note that the results in (1) and (2) remain unchanged if we include an electric 
background potential V (e.g., V bounded and measurable), while supersymmetry 
is destroyed upon addition of an electric potential. 

Under suitable assumptions, one can show ([B-1,2]) that the operators H±(Xa) 
have a common gap: in the simplest case, we assume B > 0. A more delicate 
result requires that B is positive in a neighborhood of dQ and that there exists 
a periodic gauge s.th. da — B. (By the divergence theorem, this gauge condition 
is equivalent to the flux condition f Q B(x)dx = 0). 

4. Periodic divergence type operators 

Here we report on some results of [HL-1,2] concerning periodic two-component 
media with high contrast. In the following we always assume m > 2. The 
operators dicussed below occur in acoustics, heat conduction, and propagation of 
electro-magnetic waves in photonic crystals. 
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The quadratic form is given by the Dirichlet integral 

t x [u]:= J(l + \ X n(x))\Vu(x)\ 2 dx, ueH 1 , A > 0; (4.1) 

this defines a monotonically increasing family of forms. The associated self- 
adjoint operators T A can formally be written as T x = — V • (1 + Axn)V. 

It is clear that only such functions u from the form domain survive taking 
A — > oo that are constant on each connected component of Q. If Q is connected, 
this implies that only such -u's will remain in the limiting form domain Qoo that 
are identically zero on Q. We have: 

4.1. Theorem. [HL-1] Suppose Q is open, periodic, and connected. Then T\ 
converges to the Dirichlet- Laplacian —A M on M = M. m \Q in strong resolvent 
sense. Furthermore, the associated density of states measures converge weakly. 

The standard Floquet-decomposition gives a more detailed picture. Let 
denote the operator —V • (1 + Axn)V, acting in L 2 (Q), with ^-periodic boundary 
conditions on dQ, for 9 G (— 7r, 7r] m ; cf. [RS-IV, Ku]. Then Theorem 1.1 can be 
applied in each fiber, i.e., for each fixed 9. As each has compact resolvent, 
we even obtain norm resolvent convergence in the fibers. The limiting operators 
have form domains consisting of those functions u G ^(Q) that are 
constant in each connected component of Q fl Q and satisfy ^-periodic boundary 
conditions. If Q is connected and intersects each face of the period cell Q, it 
follows that such u must vanish on Q fl Q, provided 9 ^ (0, . . . ,0). Pursuing this 
observation in detail, one obtains 

4.2. Theorem. ([HL-1]) Suppose Q C W 71 is open, periodic, connected and con- 
tains dQ. We also assume that M* (defined in Section 2) is not of measure zero. 
Let < Si < 5 2 < • • • denote the eigenvalues of —A Mn g. We then have: 

(i) The first spectral band ofT\ converges to the interval [0, Si]. 

(ii) For A large, T\ has a spectral gap following the first band. 

(Hi) The i.d.s. measure ofT\ concentrates at the set {5k~, k G N} ; as A — > oo. 

It follows from (i) that there is no norm resolvent convergence of the family 
(Tx) (The only possible candidate for the limit is —Am which has no spectrum 
below Si > 0. But norm resolvent convergence implies spectral convergence.) 

For smoothly bounded M a much more detailed picture of the process of spec- 
tral concentration is given in [F]. 

4.3. Remark. More general results can be found in [HL]. First of all, without ad- 
ditional effort one can handle elliptic operators of the type —V • (sl(x) + Xb(x)) V, 
under suitable assumptions on the (symmetric) matrices a(x) and h(x), like 
aij,bij G Loo(R m ), sl(x) uniformly positive definite, h(x) = on M, and pos- 
itive definite for x G Vt. 

Second, [HL] also gives criteria for the existence of higher gaps: if the eigenspace 
of the fc-th Dirichlet eigenvalue 5^ of — A MCi q contains an eigenfunction <i fc s.th. 
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/ d k 7^ then a(T\) will have a gap for A large that converges to (5 k ,u k ), for 
some Uk > 5k- 

4.4. Remark. There are natural and simple situations where one is led to consider 
period cells that are more general than (0, l] m . In fact, all of the above works if 
T is a discrete sublattice of R m that spans R m and if Q is a contractible subset of 
R m with piecewise smooth boundary tessalating R m via the translations 7 G T. 
Here we think, in particular, of the class of results that depend on the assumption 
that Q contains the boundary of the period cell. 

4.5. Remark. In R m , m > 3, it is possible to have both Q and M connected. A 
sponge is a typical example (albeit, in reality, not periodic but random!). Here 
our analysis would yield that, for high contrast, the i.d.s. is approximately given 
by the i.d.s. of T^. The latter can be further analyzed by Floquet-decomposition 
if we assume that Q intersects each face of Q (think of the case where M is 
connected from one cell to the neighboring cells through thin filaments only). 

4.6. Remark. It is natural to ask what happens for a decreasing sequence of 
operators like 

Sx = -V • — V, < A T 00. (4.2) 

1 + Axm 

The associated quadratic forms s A are defined by S\[u] — J(l + Axm) _1 | Vu\ 2 dx 
on ?i 1 (]R m ). Evidently, the largest form that is smaller than S\, for all A > 0, has 
domain 7i 1 (IR m ) and takes the values j n \Vu\ 2 dx. Let us call this form s^. It is 
easy to see that Soo is closable. Let denote the self-adjoint operator associated 
with the closure of Soo. It follows from Simon's theorem [RS-I; Thm. S. 16] that 
the operators S\ converge to in strong resolvent sense. 

Let us assume, for the moment, that any u G 7i 1 (fi) can be extended to a 
function in TiMW 1 ) and that C™(M°) is dense in L 2 (M). Then it is not difficult 
to show that the closure of the form has domain 7i l (Vt) © L 2 {M), and takes 
the values Soo[u] = f n \Vu(x)\ 2 dx. Under these assumptions is the direct sum 
of the Neumann Laplacian on Q and the zero operator on L^M). In particular, 
the i.d.s. of would be infinite on the positive real line, and it seems rather 
hopeless to extract any useful spectral information on S\. (Of course, the roles 
of M and VL are interchangeable here if the boundary of M is smooth enough.) 

4.7. Remark. In the theory of photonic crystals, one is led to study the "high- 
contrast /thin- wall limit" of the acoustic operator in R 2 . In typical applications, 
the optically dense medium occupies a connected, periodic set M consisting of 
thin walls, cf. [FK-1,2]. In our notation, we would then consider the limit of 
T\ = — V(l + Axq)V, as A — > 00. Note that now fl is not connected. Our results 
yield a strong resolvent limit and a Floquet-decomposition analogous to what we 
have discussed above. The limiting operators T^} in the fibers now have a form 
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domain consisting of functions that are constant on f2 H Q, where the constant 
value will be non-zero, in general. Therefore, it is much harder to determine the 
spectrum of the than in the case of Eqn. (4.1). If f2 n Q exhausts "most of" 
Q (thin wall asymptotics), then the analysis can proceed again [FK-1,2], [KK]. 

4.8. Remark. Random Media and the convergence of the i.d.s. for the operators 
of Eqn. (4.1) with a random geometry have been studied by K. Lienau [L]. Since 
we need a connected Q for our approach to work well, the ergodic process that 
generates the random set M has to respect this property. Here various methods 
from stochastic geometry ([SKM]) can be employed. 

5. Periodic Laplace-Beltrami operators 

In this section we report on some related results on periodic Laplace-Beltrami 
operators, cf. [P-1,2]. Let M. be a non-compact Riemannian manifold of dimen- 
sion d > 2. We call M. a periodic or covering manifold iff Y = IT 1 acts properly 
discontinuously and isometrically on M.. Under these conditions, the quotient 
M./Y is again a Riemannian manifold which we assume to be compact. 

We call Q a period cell or fundamental domain of M. iff Q is open, Q does 
not intersect any other translate 7 + Q except for 7 = 0, and the union of all 
translates of Q covers the whole manifold M. Note that we assume here that Q 
is an open set. One can apply Floquet decomposition in the same way as before; 
here, -^-periodicity means that u{j + x) — e^u^x) for all x G Q and 7 G T such 
that 7 + (plus an analogous condition for the first derivatives). 

Denote by (^) the metric tensor of M. in some chart. We consider conformal 
deformations of the given periodic metric, i.e., we set 

9ij(x;\) = Ql(x)gij(x). 

Here, Q\ is a strictly positive, periodic, smooth function converging pointwise to 
the indicator function of a closed periodic set M as A — > 00. In particular, we 
assume that 

(1) M := M fl Q is compactly contained in Q with smooth dM; 

(2) g x (x) = 1 for all x G M; 

(3) g x (x) = {, for all x G Q with dist(x, M) > X~ m . 

To simplify the notation, we restrict ourselves to the conformally flat case, i.e., 
we assume that M. = R m with the Euclidean metric tensor gij(x) = 5ij, defined 
on W 71 . We assume that the boundary dM has a tubular neighborhood where 
we can introduce local coordinates (r, y) given by the distance r from dM and 
y G dM. For technical reasons we need the additional assumption that these 
coordinates can be extended to Q fl Q; e.g., think of a ball M centered in the 
period cell Q = (0, l) m . The Laplace-Beltrami-operator B\ corresponding to the 
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deformed metric is given by 

V 9\ 

(using the summation convention), where g\ = det gij(-,\) = p 2 x m denotes the 
square of the volume density of the deformed metric tensor. This operator is 
defined as a non-negative, self-adjoint operator via the quadratic form 

b x [u]:= J g il {-,\)d i ud j u^dx = J \Vu\ 2 g™~ 2 dx, ueH 1 ^" 1 ) (5.1) 

in the Hilbert space ?i\ := L 2 (M m ) with inner product j uv^fg~\dx. Note that, in 
general, both the Dirichlet integral and the volume form depend on the parameter 
A. The fundamental difference to the case of divergence type operators lies in 
the dependence of the inner product on A. In particular, for m > 3, there is no 
monotonicity. 

Floquet-decomposition implies that it is enough to study ^-periodic b.v.p.'s 
on Q. Furthermore, Dirichlet-Neumann-bracketing yields enclosures of bands by 
Neumann and Dirichlet eigenvalues E^ D (X) (ordered by min-max and repeated 
according to multiplicities). Similarly, denote by the Neumann eigenvalues 
ofM = MnQ. 

Since the norm of the Hilbert space depends on A, we cannot apply Theo- 
rem 1.1. Instead, the convergence of the eigenvalues on Q to the Neumann eigen- 
values on M is obtained directly through an application of the min-max-principle 
and a careful comparison of the corresponding Raleigh quotients. 

In particular, we show that E% is (approximately) a lower bound for E^ (A) 
(using the restriction u\m as transition operator from H}{Q) to 7Y 1 (M )). In 
the same way, we prove that E^ is also an approximate upper bound for Ej?(\) 
provided m > 3 (using an extension operator from 7i 1 (Mo) to T~io(Q), the Ti 1 - 
closure of C%°(Q), as transition operator). Therefore, we have 



5.1. Theorem. ([P-1,2]) Let m > 3. 
of the Neumann Laplacian on M , as A 
A — > oo. 



Then o~(B\) converges to the spectrum 
— > oo. In particular, gaps open up, as 



5.2. Remark. The case m = 2 is special. Here, only the inner product depends 
on A. Inverting the role of the norm and the quadratic form yields a limiting 
operators in each fiber which, however, depend on 9. In the special case of a 
2-dimensional cylinder one can nevertheless prove the existence of spectral gaps 
by direct calculations (separation of variables), cf. [P-1,2]. 

5.3. Remark. In contrast to the Euclidean case, we may choose the dimension 
d of the manifold Ai independently of the dimension of periodicity m: think of 
the surface of a cylinder (d — 2, m — 1) or the surface of a periodic jungle gym 
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(d — 2, m — 3). Furthermore, one can allow any Abelian finitely generated group 
T of infinite order, i.e., products of Z and cyclic groups. 

5.4. Remark. In some sense, the case of Laplace-Beltrami operators is the oppo- 
site of the case of divergence type operators: while the divergence type operator 
needs the factor (1 + Xxn) to be large on Q to produce spectral gaps, the Laplace- 
Beltrami operator needs the factor Q™~ 2 to be small on fi. 

5.5. Remark. Davies and Harrell [DH] proved the existence of at least one gap 
in the periodic conformally flat case (transforming the conformal Laplacian into 
a corresponding Schrodinger operator). This result is a special case of our result. 
Using similar methods as [DH], Green [G] obtained examples with a finite number 
of gaps in the 2-dimensional conformally flat case. Recently, Yoshitomi [Y] proved 
the existence of spectral gaps for the (Dirichlet) Laplacian on periodically curved 
quantum wave guides. In [DH, G, Y], the existence of gaps is established by 
analyzing a one-dimensional problem. Here, in contrast, we directly study the 
multi-dimensional problem. 

5.6. Remark. More "geometric" constructions of periodic manifolds with spectral 
gaps can be found in [P-1,2], think for example of an infinite number of copies 
of a compact manifold, joined by thin cylinders. If these cylinders are small 
or thin enough, one again obtains spectral gaps of the corresponding Laplace- 
Beltrami operator. These results are closely related to earlier works of Chavel 
and Feldman [CF] and Anne [A] . 

6. Maxwell operators 
Another important elliptic operator is the Maxwell operator given as 

M = —curl— 3— curl, (6.1) 

e(x) 

acting on solenoi'dal vector fields in IR 2 or in IR 3 , where e(x) denotes the dielec- 
tricity of the medium. 

The existence of gaps for periodic structures is of fundamental importance for 
photonic crystals. In many cases, one is satisfied to find intervals of low spectral 
density (quasi-gaps). This class of problems has been studied in depth in a series 
of papers by A. Figotin, P. Kuchment and others (cf., e.g., [FK-1,2], [KK]) where 
a complete analysis for high-contrast and/or thin- wall asymptotics is given. 
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